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2004-2005  ANNUAL  PROGRESS  REPORT 


This  report  presents  the  specific  aims  and  accomplishments  of  our  breast  cancer  research  project 
during  the  first  year  of  funding  sponsored  by  the  US  Department  of  the  Army. 


Introduction 

The  goal  of  this  research  project  is  to  apply  the  multiple  monitoring  techniques,  i.e.  Near  infrared 
spectroscopy  (NIRS),  FOXY  oxygen  sensor  and  ,9F  MR  EP  imaging  of  Hexafluorobenzene  (HFB), 
to  prove  the  following  hypotheses:  combination  of  hyperbaric  oxygen  (HBO)  intervention  can 
significantly  improve  breast  tumor  oxygenation,  and  that  tumor  oxygenation  remains  elevated  for  a 
substantial  period  of  time  even  after  HBO  exposure,  which  may  be  a  novel  approach  to  enhance 
radiosensitivity.  If  our  hypotheses  are  proven  to  be  true,  this  study  will  lead  to  an  optimal 
intervention  plan  to  improve  tumor  oxygenation  and  to  determine  an  optimal  time  interval  after 
HBO  decompression  for  radiotherapy.  Such  a  novel  approach  will  largely  enhance  the  efficiency  of 
non-surgical  therapies  for  breast  tumor  treatment  and  provide  a  novel  prognostic  tool  for  clinical 
practice.  This  study  will  also  provide  a  better  understanding  of  tumor  vasculature  and  tissue  oxygen 
dynamics  and  spatial  heterogeneity  under  HBO  exposure. 

The  project  has  three  specific  aims: 

Aim  1:  to  determine  the  absolute  values  of  oxygenated  hemoglobin  concentration,  [Hb02],  and 
hemoglobin  oxygen  saturation,  SOi,  in  solid  breast  tumors  from  the  NIRS  measurements. 

Aim  2:  to  investigate  vascular  oxygenation  and  tissue  oxygen  tension  of  breast  tumors  under 
continuous  normobaric  and  hyperbaric  oxygen  exposures  with  several  gas  interventions,  using 
both  a  single-channel  NIRS  system  and  3-channel  FOXY  p02  system  simultaneously. 

Aim  3:  to  investigate  global  and  local  dynamics  of  tumor  vascular  [Hb02]  and  tissue  DO2  of 
breast  tumors  immediately  after  HBO  exposure  by  using  both  three-channel  NIRS  and  1  r  MR 
EP  imaging  simultaneously. 

Body  of  the  report 

As  mentioned  above,  the  purpose  of  this  project  is  to  investigate  vascular  and  tissue  oxygen 
dynamics  in  breast  tumors  by  using  multiple  techniques:  NIRS,  FOXY  and  MRI.  In  the  first  year, 
my  focus  is  to  develop  algorithm  to  calculate  absolute  value  of  hemoglobin  concentration  (month  1- 
8),  to  obtain  the  appropriate  knowledge  and  skills  for  conducting  the  project  (month  9-14).  1  have 
accomplished  not  only  those  two  tasks  (month  1-14)  but  also  a  part  of  task  3  (month  14-30). 

Task  1:  to  develop  a  numerical  algorithm  so  as  to  determine  the  absolute  value  of  HhO?  and 
SO 2  in  solid  breast  tumors  from  near  infrared  spectroscopic  (NIRS)  measurements  by  using 
the  non-linear,  least-squared  optimization  apyroach.  (months  1-8) 

Since  task  1  has  been  accomplished  by  other  members  of  our  group,  the  PI  has  switched  to 
another  task  of  equal  importance. 
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New  task  1:  develop  a  mathematical  model ,  which  relates  the  blood  flow  and  metabolic  rate  of 
oxygen  in  breast  tumors  to  hemoglobin  concentration  so  that  NIRS  can  be  used  to  estimate 
changes  in  tumor  blood  flow  (TBF)  and  tumor  metabolic  rate  of  oxygen  (TMRO?)  non- 
invasivelv  during  hvveroxic  gas  interventions  in  rat  tumors.  Such  a  model  will  shed  light  on 
tumor  vascular  adaptation  to  therapeutic  intervention,  giving  us  a  better  understanding  of 
physiological  responses  of  tumors  to  intervention.  This  task  has  been  accomplished,  as  reported 
below: 

In  order  to  find  the  relationship  between  the  hemoglobin  concentration  and  metabolic  rate  of 
oxygen  during  intervention,  there  is  a  mathematical  equation,  called  as  the  ratio  method  which 
states  that  metabolic  rate  of  oxygen  could  be  derived  from  blood  flow  and  hemoglobin 
concentration  (Jones  et  al  2001).  NIRS  measurements  have  enabled  us  to  quantify  relative 
changes  in  hemoglobin  concentration  during  inhalation  but  have  not  been  used  to  measure  blood 
flow  changes.  Thus,  it  is  necessary  to  explore  a  mathematical  model  that  relates  changes  in  blood 
flow  to  changes  in  blood  volume  (i.e.,  total  hemoglobin  concentration)  in  order  to  estimate  flow 
changes  based  on  the  NIRS  data.  Boas  et  al  have  utilized  the  Windkessel  model  to  fit  a  flow- 
volume  relationship  and  estimated  the  relative  changes  of  cerebral  blood  flow  during  brain 
activation  (Boas  et  al  2003),  with  a  relatively  good  fit  to  the  NIRS  experimental  data.  Following 
the  same  strategy,  here,  oxygen/carbogen  induced  tumor  blood  flow  changes  were  estimated  by 
fitting  the  Windkessel  model  (Mandeville  et  al  1999)  to  the  measured  A[HbT],  with  the 
modification  of  physiological  meaning  of  tumor  vasculature. 

Step  1:  The  Windkessel  model  used  for  estimating  blood  flow  from  [HhTl 

The  Windkessel  model  relates  blood  flow  to  blood  volume.  With  a  few  assumptions,  we  can  thus 
estimate  the  blood  flow  response  to  intervention  from  the  measured  A[HbT]  since  A[HbT]  is 
directly  proportional  to  changes  in  blood  volume.  The  Windkessel  model  is  based  on  the 
conservation  of  mass  to  associate  changes  in  blood  volume  with  changes  in  blood  flow  through 
the  regional  arterial,  capillary  and  venous  compartments,  with  capillaries  and  veins  lumped 
together  and  named  a  Windkessel  compartment.  The  model  proposes  that  the  flow  into  the 
Windkessel  compartment  is  largely  determined  by  the  vasomotor  control  of  arterioles,  and  that 
the  capillaries  and  veins  passively  respond  to  arterial  pressure  changes  (Mandeville  et  al  1999). 
Briefly  described  here,  the  model  has  the  following  definitions  and  relationships: 

1)  In  analogy  to  Ohm’s  law,  we  have  P(t)=F(t)R(t), 

where  F  is  the  blood  flow,  P  is  the  blood  pressure,  and  R  is  the  vascular  resistance; 

2)  The  Windkessel  model  assumes  Vw(t)=APw(t)'^, 

where  Vw  and  Pw  are  the  volume  and  pressure  of  the  Windkessel  compartment, 

respectively. 

3)  The  Windkessel  model  also  assumes  Rw(t)/Rw(0)=(V(0)/Vw(t))a, 
where  Rw  is  the  resistance  in  the  Windkessel  compartment. 

4)  The  arterial  volume  changes  are  related  to  the  arterial  resistance  changes. 
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In  the  relationships  given  above,  (3  represents  the  vascular  compliance,  a= 2  indicates  laminar 
flow  with  the  vessel,  A  is  a  constant  and  equal  to  Vw ( 0)/(F( O)Rw(O))  1 11 .  With  these  definitions 
and  initial  conditions,  the  Windkessel  model  then  arrives  at  the  coupled  differential  equations, 
(1)  and  (2),  for  flow  and  volume  changes  due  to  arterial  resistance  changes,  i.e. 


dVJQ 

81 


VJ£T^_ 

0)K(0)“ 


(1) 


P-Pjl)  P-Vs,(lflA,> 
PAD  ~  PA') 


(2) 


In  this  study,  we  assume  a  bi-phasic  model  for  the  temporal  response  of  arterial  resistance  during 
hyperoxic  gas  intervention  according  to  the  bi-phasic  behavior  of  hemodynamic  responses  to 
oxygen  or  carbogen  intervention,  i.e. 

Ra  =  /?^(0)exp(-^)  +  7^s(0)exp(-^)  (  3  ) 

T.f  Ts 


where  Ra/ and  RA.,  are  the  amplitudes  of  fast  and  slow  arterial  resistance  components,  T/  and  rv  are 
the  time  constants  of  fast  and  slow  component,  respectively,  and  to  is  the  time  when  the  rats  were 
exposed  to  the  hyperoxic  gas  intervention.  The  numerical  solutions  of  these  differential 
equations,  (1)  and  (2),  then  provide  the  relationship  between  blood  flow  and  blood  volume.  The 
units  are  normalized  as  F,„(0)=F„,„(0)=1,  Vu(0)=F,„(0)t,  and  Wi(0)+i?M(0)=l . 


In  order  to  estimate  the  blood  flow  from  blood  volume  changes,  we  must  assume  a  constant 
value  of  hematocrit  during  the  intervention  so  that  we  will  have  the  relationship  of 
A[HbT]/[HbT]o=ABV/BVo.  Based  on  the  Windkessel  model,  four  parameters  for  the  arterial 
resistance,  i.e.,  Ra/ 0),  Ra/0),  x/and  rv,  need  to  be  quantified  besides  the  Windkessel  vascular 
reserve  p  and  the  Windkessel  transit  time  x.  The  model  parameters  were  determined  or  optimized 
in  a  nonlinear  fit  to  the  experimentally  measured  A[HbT],  Given  the  optimal  parameters,  the 
Windkessel  model  permits  to  estimate  ATBF/TBFo,  which  is  the  function  of  arterial  resistance 
(Ra),  p  and  x. 


Step  2:  Ratio  method  to  estimate  A  TMRO? 

Ratio  method  has  been  used  to  calculate  the  change  in  cerebral  metabolic  rate  of  oxygen  (Jones 
et  al  2001).  Here  we  will  use  the  same  method,  with  the  C  (cerebral)  replaced  by  T  (tumor)  in 
equation  4: 


0+ 


atmro2 

TMRO , 
*0 


)  =  (!  + 


A  TBF,  „  Mib 

- Ml  +  y - 

TBF0  Hb„ 


)/(l  +  Z7- 


A HbT 
HbT„ 


(4) 


where  y  and  jj  represent  the  relative  changes  of  deoxy-hemoglobin  and  total  hemoglobin 
concentration  in  all  vascular  compartments,  respectively.  In  our  study,  A[HbT]  and  A[Hb]  are 
measured  by  NIRS,  while  ATBF  can  be  estimated  from  the  Windkessel  model.  Then,  eq.  (4)  will 
further  allow  us  to  compute  an  estimate  for  changes  in  TMRO2. 
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Step  3:  Estimation  of  blood  flow  and  ATMRO?  of  rat  breast  tumors  under  intervention 

There  were  totally  1 1  rats  under  investigation.  Group  1  (n=6)  was  under  air-oxygen  inhalation, 
shile  group  2  (n=5)  had  air-carbogen  intervention.  Figure  1  shows  an  example  of  time  course  of 
tumor  A[HbT]  and  A[Hb02].  Raj  arid  RAs,  the  amplitudes  of  fast  and  slow  arterial  resistance 
change  in  response  to  interventions,  are  equal  to  0.19  and  0.6  respectively.  The  two  thick  curves 
were  obtained  by  averaging  six  sets  of  measurements  from  the  rat  tumors  in  group  1  when  the 
rats  were  subject  to  oxygen  intervention.  Both  A[HbT]  and  A[Hb02]  have  stable  baselines  when 
the  rats  were  breathing  air.  When  the  gas  was  switched  to  oxygen,  both  parameters  increased 
consistently  with  the  typical  bi-phasic  pattern.  ATBF  was  estimated  from  A[HbT],  based  on  the 
Windkessel  model  as  shown  in  eqs.  (1)  to  (3);  then,  ATMRCF  was  obtained  using  eq.  (4).  Both 
ATBF  and  ATMR02  have  exhibited  significant  increases  with  a  bi-phasic  pattern  too. 


0  10  20  30  40 

Time  (min) 


Figure  1.  Time  course  of  averaged  values  of  A[HbT]  (light  thick  curve),  A[Hb02]  (dark 
thick  curve),  with  fitted  ATBF  (thin  dotted  curve)  and  ATMRO?  (thin  dark  curve),  taken 
from  six  tumors  with  gas  intervention  of  air-oxygen.  The  units  of  A[HbT]  and  A[Hb02] 
are  mM/DPF,  and  the  fitted  ATBF  and  ATMR02  have  arbitrary  units. 

However,  the  tumors  have  various  responses  of  A[HbT],  ATBF,  and  ATMR02  to  carbogen 
inhalation  (Fig  2  to  4),  while  tumor  A[Hb02]  always  gave  consistent  increases  with  given 
carbogen  intervention  (Fig  2  to  4).  Figure  2  shows  the  time  course  of  averaged  A[Hb02], 
A[HbT],  computed  ATBF  and  ATMRO2  from  three  tumors  in  group  2  (tumor  #1~3)  when  the 
rats  were  subject  to  the  gas  intervention.  In  this  case,  there  is  a  consistent  increase  of  A[Hb02] 
and  A[HbT],  Estimated  from  A[HbT]  using  the  Windkessel  model,  Raj  and  RA,X  are  equal  to  0.2 
and  0.6  respectively,  ATBF  was  elevated,  and  so  was  the  corresponding  ATMR02.  However, 
A[HbT]  taken  from  Tumor  #4  from  group  2  did  not  show  much  change,  and  neither  did  ATBF 
and  ATMRO2  from  this  tumor  in  response  to  carbogen  inhalation  (Figure  3),  while  the  fitted  RAf 
and  Ra,s  are  equal  to  0.05  and  0.05  respectively.  Furthermore,  Tumor  #5  in  group  2  shows  a 
decreased  A[HbT]  induced  by  carbogen  intervention,  with  reduced  ATBF  and  ATMRO2 
accordingly  (Figure  4),  and  the  fitted  RAjand  RA  s  are  equal  to  0.45  and  0.5  respectively. 
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Figure  2.  Time  course  of  average  values  of  A[HbT]  (light  thick  curve),  A[Hb02]  (dark 
thick  curve),  with  fitted  ATBF  (thin  dotted  curve)  and  ATMR02  (thin  dark  curve),  from 
three  tumors  with  gas  interventions  of  air-carbogen.  The  units  of  A[HbT]  and  A[Hb02] 
are  mM/DPF,  and  the  fitted  ATBF  and  ATMR02  have  arbitrary  units. 
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Figure  3.  Time  course  of  A[HbT],  A[Hb02],  computed  ATBF  and  ATMR02  from 
Tumor  4  of  group  2  with  air-carbogen  intervention. 
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Figure  4.  Time  course  of  A[HbT],  A[Hb02],  computed  ATBF  and  ATMRCL  from 
Tumor  5  of  group  2  with  air-carbogen  intervention. 


Task  2:  to  obtain  the  knowledge  and  appropriate  training  for  handline  laboratory  animals  and 
performing  the  measurements,  (months  9-14) 

a.  To  have  training  in  handling  laboratory  animals. 

b.  To  study  the  fundamentals  and  operation  of  the  1 -channel  and  3-channel  TMIRS  system. 

c.  To  study  the  fundamentals  and  operation  of  FOXY  fiber  optic  oxygen  sensor. 

d.  To  study  the  fundamentals  and  operation  of  19F  MR  EP1  imaging. 

e.  To  study  the  fundamentals  and  operation  of  hyperbaric  chamber. 

The  PI  has  obtained  the  knowledge  and  training  for  handling  laboratory  animals  and  performing  the 
measurements.  The  PI  has  also  studied  the  fundamentals  and  operations  of  multi-channel  FOXY 
oxygen  sensor,  multi-channel  broadband  diffuse  NIR  spectroscopy,  and  l9F  EP  imaging  and 
hyperbaric  chamber. 

Task  3:  to  study  the  influence  of  four  different  sas  interventions  on  tumor  oxygenation  using 
both  single-channel  NIRS  and  FOXY  system,  on  various  rat  breast  tumor  size,  (months  14-30) 

The  PI  has  made  tremendous  efforts  to  have  the  fiber  probes  ready  for  both  the  NIR  system  and 
FOXY  oxygen  sensor  as  well  as  to  integrate  the  probes  with  the  hyperbaric  chamber.  We 
implemented  single-channel  NIRS  with  broadband  diffuse  spectroscopy  (BDS)  because  BDS 
provides  a  broad  spectrum  from  400  to  900  nm,  enabling  us  to  calculate  the  absolute  values  of 
hemoglobin  concentration  with  the  algorithm  that  we  are  currently  developing.  Being  portable,  BDS 
consists  of  a  tungsten  halogen  light  source  (HL-2000HP,  ocean  optics,  FL),  a  hand-held 
spectrometer  (USB2000,  Ocean  optics,  FL),  and  two  optical  fiber  bundles  with  respective  2.6  mm 
and  1  mm  bundle  diameter.  The  current  design  of  the  fibers  provides  us  with  better  optical 
transmittance  without  sacrificing  their  compatibility  of  light  source  and  spectrometer. 
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Figure  5.  Picture  of  hyperbaric  chamber  (left  one), 
integrating  with  broadband  light  diffuse  spectroscopy 
source  and  detector  probes  (white  ones)  and  FOXY  oxygen 
sensor  probes  (blue  ones),  FOXY  oxygen  sensor  (right  one) 
and  BDS  (middle  ones) 


As  shown  in  Figure  5,  the  hyperbaric  chamber  is  made  of  transparent  acrylic,  which  allows  direct 
visualization  during  experiments.  The  chamber  is  cylinder-shaped  with  an  inner  diameter  of  8  inch 
and  a  length  of  24  inch.  One  end  is  attached  with  a  removable  door,  and  the  other  end  is  fixed,  on 
which  a  pressure  meter  and  gas  inlet  and  outlet  valves  are  integrated  to  monitor  and  control  the 
pressure  and  gas  flow  rate  in  the  chamber.  Optical  fiber  probes  of  BDS  and  FOXY  were  assembled 
into  the  chamber,  while  the  instruments  themselves  were  placed  outside  the  chamber  to  avoid  any 
possible  fire  hazards  (Figure  5). 

Because  safety  is  a  very  important  consideration  for  hyperbaric  chamber  (Jain  1999),  all  facilities 
required  in  the  chamber  are  specially  designed  to  avoid  any  possibility  of  fire  hazard.  The  optical 
probe  ends  were  made  of  plastic  instead  of  metal  to  avoid  static  electricity,  which  may  be  produced 
by  the  friction  between  the  probe  and  probe  holder.  The  homemade  optical  breadboard  inside  the 
chamber  is  made  of  acrylic  plastic,  with  two  fixing  tubes  attached  on  both  sides  to  assure  the 
stability  of  the  board  sitting  in  the  cylindrical  chamber.  All  metal  probe  holders  were  washed 
thoroughly  to  remove  all  lubricate  oil. 


Figure  6.  Schematic  diagram  of  experimental  setup  for  simultaneous  monitoring  of 
tumor  vascular  and  tissue  oxygenation  when  the  rats  are  under  hyperbaric  gas 


A  Fisher  344  rat  was  placed  in  the  hyperbaric  chamber  before  the  experiment,  as  shown  in  figure  6. 
The 


probe  ends  of  source  and  detector  bundles  of  BDS  were  placed  on  the  opposite  side  of  the  tumor  in 
transmission  geometry.  Three  needle  probes  of  the  FOXY  oxygen  sensor  were  inserted  into 
different  regions  of  the  tumor.  The  rat  was  exposed  to  varying  gas  environment  according  to  the 
following  challenge  paradigms: 

Air(15min)  Oxygen  (1  ATA,  15  min)  ->  Oxygen  (2  ATA,  30  min)  ->  Air  (30  min) 
where  ATA  stands  for  1  atmosphere  absolute. 
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Figure  7  Temporal  profiles  of 
A[Hb02]  and  p02  in  response  to 
respiratory  challenge  for  a  rat  breast 
tumor,  measured  simultaneously  by 
both  BDS  and  FOXY  oxygen  sensor. 
P 0;n ^ . e r  p O siavei^pf^2  siave2  represent 
p02  readings  from  three  different 
regions  of  the  tumor  (region  1 , 2,  3). 


Figure  7  shows  A[HbC>2]  and  pC>2  responses  to  respiratory  challenge  simultaneously  measured  by 
BDS  and  FOXY,  respectively.  Baseline  measurements  (breathing  air)  were  generally  stable  in 
A[Hb02],  and  altering  the  inhaled  gas  to  oxygen  in  1  atmosphere  absolute  (ATA)  induced  a  rapid 
and  significant  increase,  followed  by  a  gradual  change  in  A[Hb02].  Changing  the  gas  to  hyperbaric 
oxygen  in  2  ATA  leads  to  additional  improvement  of  A[Hb02].  Upon  return  to  air  (baseline), 
A[Hb02]  dropped  quickly  first  and  more  gradually  later.  Comparatively,  not  all  regional  p02  show 
stable  baseline:  p02  in  region  1  shows  the  stable  value  in  the  first  10  min,  but  it  has  a  gradual 
decrease  in  the  following  5  min,  while  p02  in  regions  2  and  3  remain  unfaltering  in  air.  When  the 
gas  was  switched  to  oxygen  (1  ATA),  pC>2  in  region  1  increase  gradually,  while  pCF  in  regions  2 
and  3  do  not  show  much  significant  improvement.  Further  changing  the  gas  to  hyperbaric  oxygen  (2 
ATA)  causes  substantial  improvement  of  pC>2  in  all  three  local  regions.  Upon  return  to  air,  pC>2  in 
regions  1  and  2  return  to  baseline  gradually,  while  that  in  region  3  drops  to  baseline  quickly. 


Key  Research  Accomplishments 

1)  The  PI  had  developed  a  mathematic  model  to  correlate  the  hemoglobin  concentration  to  blood 
flow  and  metabolic  rate  of  oxygen  in  tumor.  We  demonstrated  that  we  are  able  to  quantify  the 
tumor  blood  flow  and  metabolic  rate  of  oxygen  based  on  the  measurement  of  hemoglobin 
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concentration  using  NIRS  during  gas  intervention. 

2)  The  PI  has  obtained  the  knowledge  and  training  for  handling  laboratory  animals  and  performing 
initial  animal  measurements  independently  using  multi-channel  FOXY  oxygen  sensors, 
broadband  diffuse  NIR  spectroscopy,  19F  EP1  imaging  and  a  hyperbaric  chamber. 

3)  The  PI  has  designed  and  implemented  fiber  probes  ready  for  both  the  NIR  system  and  FOXY 
p02  system,  and  integrated  the  probes  with  the  hyperbaric  chamber.  She  had  gotten  preliminary 
data  of  the  simultaneous  measurement  under  hyperbaric  oxygenation,  which  is  supposed  to 
accomplish  in  the  second  year. 

Reportable  Outcomes: 

Peer-reviewed  Journal: 

1)  Mengna  Xia,  Vikram  Kodibagkar,  Hanli  Liu  and  Ralph  Mason,  “Tumor  oxygen  dynamic 
measured  simultaneously  by  near  infrared  spectroscopy  and  19F  magnetic  resonance  imaging  in 
rats”,  under  2nd  review  of  Physics  in  Medicine  and  Biology. 

2)  Jae  Kim,  Mengna  Xia,  Hanli  Liu,  “Extinction  coefficients  of  hemoglobin  for  near-infrared 
spectroscopy  of  tissue”,  IEEE  Engineering  in  Medicine  and  Biology  magazine,  24:  118-121 
(2005) 

3)  Yueqing  Gu,  Wei  R.  Chen,  Mengna  Xia,  Sang  W.  Jeong,  and  Hanli  Liu,  “  Effect  Of 
photothermal  therapy  On  breast  tumor  vascular  contents:  non-invasive  monitoring  by  near 
infrared  spectroscopy”,  Photochemistry  and  Photobiology,  in  press,  2005 

Presentation  and  Proceeding  Papers: 

1)  Mengna  Xia,  Vikram  Kodibagkar,  Ralph  Mason,  Benjamin  Levine,  Hanli  Liu,  “Tumor  vascular 
and  tissue  oxygenation  dynamics  under  normobaric  and  hyperbaric  oxygen  interventions”, 
presented  at  the  fourth  Era  of  Hope  meeting  for  the  Department  of  Defense  (DOD)  Breast 
Cancer  Research  Program  (BCRP)  held  on  June  8-11, 2005  in  Philadelphia,  Pennsylvania 

2)  Mengna  Xia,  Ralph  Mason,  Hanli  Liu,  “A  model  of  hemodynamic  responses  of  rat  tumors  to 
hyperoxic  gas  challenge”,  Proc.  SPIE-  Optical  Tomography  and  Spectroscopy  of  Tissue  VI, 
5693:301-307(2005) 

3)  Mengna  Xia,  Vikram  Kodibagkar,  Yueqing  Gu,  Anca  Constantinescu,  Ralph  Mason,  Hanli  Liu, 
“Tumor  oxygen  dynamics  measured  simultaneously  by  near  infrared  spectroscopy  and  19F  MR 
EPI  imaging”,  in  Biomedical  Topical  Meetings  on  CD-ROM  (The  Optical  Society  of  America, 
Washington,  DC,  2004),  ThF33 

4)  Yueqing  Gu,  Vikram  Kodibagkar,  Mengna  Xia,  Anca  Constantinescu,  Ralph  P.  Mason,  Hanli 
Liu,  “Breast  tumor  vascular  oxygenation  and  blood  volume  assessed  by  near-infrared 
spectroscopy  and  magnetic  resonance”,  Proc.  SPIE-  Optics  in  Health  Care  and  Biomedical 
Optics:  Diagnostics  and  Treatment  II,  5630:  918-925  (2004) 

Conclusions: 

Up  to  date,  we  could  draw  the  following  conclusions  from  the  work  that  we  have  conducted  in 

Year  1: 

1)  We  estimated  the  hyperoxic  gas  induced  tumor  blood  flow  changes  by  fitting  the  model  to 
the  measured  A[HbT].  The  bi-phasic  feature  of  the  fitted  blood  flow  is  consistent  with  the 
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hypothesis  that  there  exist  two  distinct  regions  within  a  tumor:  1)  a  well  perfused  region  with 
a  regular  blood  flow,  and  2)  a  poorly  perfused  region  with  slow  blood  flow  (***  references 
of  my  Applied  Optics  paper  and  Jae’s  Optics  Express  paper).  Mathematically,  the  bi-phasic 
feature  of  blood  volume  (i.e.,  total  hemoglobin  concentration)  change  during  hyperoxic  gas 
intervention  provides  a  connection  between  the  computational  model  and  our  experimental 
tumor  data.  We  believe  that  the  hemodynamic  parameters,  such  as  A[Hb02],  A[Hb]  and 
A[HbT]  directly  measured  from  N1RS,  ATBF  and  ATMR02  indirectly  estimated  from  NIRS, 
can  provide  valuable  insight  into  the  tumor  compartment  of  such  a  complex  system  and  allow 
us  to  explore  dynamic  signatures. of  tumours.  All  of  these  can,  in  turn,  enhance/assist  human 
cancer  diagnosis  and  prognosis.  Overall,  NIRS  is  a  noninvasive  tool  to  monitor  tumor 
oxygenation  and  may  also  serve  as  a  mediator  to  estimate  tumor  perfusion  and  oxygen 
consumption  rate  with  response  to  intervention  in  tumor  tissue. 

2)  Normobaric  hyperoxic  gas  (oxygen  and  carbogen)  can  effectively  manipulate  the  vascular 
oxygen  level  in  breast  tumors.  However,  different  regions  in  a  tumor  show  heterogenous 
responses  to  normobaric  oxygen  inhalation.  Hyperbaric  oxygen  is  shown  to  improve  both  the 
tumor  vascular  and  tissue  oxygen  level  more  effectively  than  normobaric  oxygen. 
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Tumor  oxygen  dynamics  by  simultaneous  NIRS  and  ,9FMR1  in  rats 


Abstract 

Simultaneous  near-infrared  spectroscopy  (NIRS)  and  magnetic  resonance  imaging  (MRI)  were  used  to 
investigate  the  correlation  between  tumour  vascular  oxygenation  and  tissue  oxygen  tension  dynamics  in 
rat  breast  13762NF  tumours  with  respect  to  hyperoxic  gas  breathing.  NIRS  directly  detected  global 
variations  in  the  oxygenated  haemoglobin  concentration  (A[Hb02j)  within  tumours  and  oxygen  tension 
(p02)  maps  were  achieved  using  l9F  MRI  of  the  reporter  molecule  hexafluorobenzene.  Multiple 
correlations  were  examined  between  rates  and  magnitudes  of  vascular  (A[Hb02])  and  tissue  (p02) 
responses.  Significant  correlations  were  found  between  response  to  oxygen  and  carbogen  breathing  using 
either  modality.  Comparison  of  results  for  the  two  methods  showed  a  correlation  between  the  vascular 
perfusion  rate  ratio  and  the  mean  p02  values  (R2>0.7).  The  initial  rates  of  increase  of  A[Hb02]  and  the 
slope  of  dynamic  p02  response,  d(p02)/dt  of  well  oxygenated  voxels  in  response  to  hyperoxic  challenge 
were  also  correlated.  These  results  demonstrate  the  feasibility  of  simultaneous  measurements  using  NIRS 
and  MRI.  As  expected,  the  rate  of  p02  response  to  oxygen  is  primarily  dependent  upon  the  well-perfused 
rather  than  poorly  perfused  vasculature. 

Keywords:  Near  infrared  spectroscopy,  Magnetic  resonance  imaging,  Breast  tumour,  Oxygenation, 
Haemoglobin,  Hexafluorobenzene. 


3 


Tumor  oxygen  dynamics  by  simultaneous  NIRS  and  19 F  MRI  in  rats 


1.  Introduction 

Tumour  oxygenation  has  been  widely  recognized  as  a  pivotal  factor  in  the  efficacy  of  radiotherapy  (Hall 
1994),  photodynamic  therapy  (Chapman  et  al  1991)  and  some  chemotherapies  (Brown  1999)  and  patient 
stratification  with  respect  to  tumour  oxygenation  status  could  be  clinically  important  (Hockel  et  al  1996, 
Fyles  et  al  1998,  Welch  et  al  2003).  It  has  been  hoped  that  modulation  of  tumour  oxygenation  could  be 
applied  to  enhance  therapeutic  efficacy.  An  attractive  intervention  is  breathing  hyperoxic  gas,  and  indeed, 
several  clinical  trials  have  examined  the  efficacy  of  normobaric  or  hyperbaric  oxygen,  to  improve 
therapeutic  outcome,  but  often  with  marginal  success  (Overgaard  and  Horsman  1996).  It  has  been 
suggested  that  outcome  might  have  been  improved,  if  responsive  tumours  could  have  been  identified  a 
priori.  Accordingly,  accurate  evaluation  of  tumour  oxygenation  in  response  to  interventions  at  various 
stages  of  growth  should  provide  a  better  understanding  of  tumour  response  to  therapy,  potentially 
allowing  therapy  to  be  tailored  to  individual  characteristics. 

Given  the  importance  of  tumour  oxygenation,  many  techniques  have  been  developed  based  on 
microelectrodes,  optical  reflectance,  Electron  Paramagnetic  Resonance  (EPR),  Magnetic  Resonance 
Imaging  (MRI)  and  nuclear  medicine  approaches,  as  reviewed  previously  (Mason  et  al  2002,  Zhao  et  al 
2004).  While  each  approach  has  unique  strengths,  some  are  highly  invasive.  Near-infrared  spectroscopy 
(NIRS)  has  been  developed  in  recent  years  as  a  promising  non-invasive  technique  to  quantify  the 
concentration  of  tissue  chromophores,  such  as  oxygenated  and  deoxygenated  haemoglobin,  water  and 
lipid  (Sevick  et  al  1991 ,  Liu  et  al  2000).  Due  to  the  deep  penetration  depth  and  biochemical  specificity  of 
NIRS,  it  has  been  widely  applied  for  quantitative  measurements  of  cerebral  oxygenation  (Delpy  and  Cope 
1997,  Yodh  and  Boas  2003)  and  blood  oxygenation  in  muscles  in  vivo  (Homma  et  al  1996).  Recently, 
NIRS  has  been  also  used  to  monitor  tumour  vascular  oxygenation  with  respect  to  interventions  (Hull  et  al 
1999,  Liu  et  al  2000,  Gu  et  al  2003).  With  notable  exceptions  (Chance  1997,  Dehghani  et  al  2004,  Shah 
et  al  2004,  Pogue  et  al  2003),  NIRS  currently  lacks  spatial  resolution,  and  thus,  the  utility  of  global 
measurements  requires  validation,  given  the  well-documented  heterogeneity  of  tumour  oxygenation.  In 
this  regard,  Conover  et  al  (2000)  compared  the  spatially  averaged  measurement  of  tumour  oxygen 
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saturation  (S02)  using  NIRS  with  the  local  S02  reading  in  individual  blood  vessels  measured  by 
cryospectrophotometry.  The  sensitivity  and  specificity  analysis  suggests  that  NIRS  may  identify  clinically 
relevant  hypoxia,  even  when  its  spatial  extent  is  below  the  resolution  limit  of  the  NIRS  technique.  We 
have  previously  investigated  correlates  between  p02  assessed  by  electrodes  (Mason  et  al  2003,  Kim  et  al 
2003b)  or  fibre-optic  probes  and  NIRS  (Gu  et  al  2003).  On  occasion,  there  was  a  good  correlation 
between  global  vascular  oxygenation  assessed  and  local  p02  at  individual  locations,  but  often,  disparate 
behaviour  was  observed.  Sequential  MRI  and  NIRS  suggested  a  better  relationship  based  on  average  p02 
from  multiple  locations  (Kim  et  al  2003b).  We  have  now  implemented  simultaneous  NIRS  and  19F  MRI 
to  examine  the  relationships  further.  Since  NIRS  is  entirely  non-invasive  it  would  provide  an  attractive 
surrogate  for  monitoring  tumour  oxygenation,  and  hence,  we  seek  correlations  with  absolute  p02 
measurements  observed  simultaneously  by  MRI. 

2.  Materials  and  methods 

Investigations  were  approved  by  the  Institutional  Animal  Care  and  Use  Committee. 

2. 1.  Animal  preparation  and  experimental  set-up 

Mammary  adenocarcinomas  13762NF  (cells  originally  obtained  from  the  Division  of  Cancer 
Therapeutics,  NCI)  were  implanted  in  skin  pedicles  (Flahn  et  al  1993)  on  the  foreback  of  ten  adult  female 
Fisher  344  rats  weighing  — 1 50  g.  When  the  tumours  reached  ~1  cm  in  diameter,  the  rats  were 
anaesthetised  with  ketamine  hydrochloride  i.p.  (100  mg/kg  body  weight,  Aveco,  Fort  Dodge,  I  A)  and 
were  maintained  under  general  gaseous  anaesthesia  (air  and  1%  isoflurane;  Baxter  International  Inc, 
Deerfield,  IL).  Tumour  hair  was  trimmed  to  give  good  optical  contact  for  NIR  light  transmission. 
Hexafluorobenzene  (HFB,  50  pi,  99.9%,  Lancaster  Co.,  Pelham,  NH)  was  administered  along  two  or 
three  tracks  in  central  and  peripheral  regions  of  the  tumours  in  a  single  plane  (transverse  to  the  rat’s 
tumour,  and  in  the  region  of  NIR  photon  pathway)  using  a  Hamilton  syringe  with  a  32  G  needle.  The 
needle  was  inserted  manually  to  penetrate  across  the  tumour  and  was  withdrawn  ~1  mm  to  reduce 
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pressure  and  3~4  pi  of  HFB  were  deposited.  The  needle  was  then  repeatedly  withdrawn  1~2  mm  and 
further  HFB  deposited  at  each  point,  as  described  in  detail  previously  (Zhao  et  al  2004). 

The  tumour  was  placed  inside  a  size-matched  Helmholtz  coil,  specially  designed  for  the  simultaneous 
MRI-NIRS  study.  The  tumour  was  inserted  between  the  two  loops  of  the  Helmholtz  coil  and  two  NIRS 
probes  were  introduced  through  the  ends  of  the  coil  along  the  coil  axis  (figure  1).  The  probes  were 
positioned  to  be  in  the  same  plane  as  the  HFB  injection.  The  rats  were  placed  in  the  magnet  on  their  side, 
and  body  temperature  was  maintained  using  a  warm  water  blanket.  A  total  of  ten  rats  were  used  in  the 
study:  seven  rats  were  subjected  to  respiratory  challenge  in  the  sequence  of  air-oxygen-air-carbogen-air, 
one  rat  breathed  air-carbogen-air-oxygen-air,  one  rat  breathed  air-carbogen-air,  and  the  last  one  breathed 
air-oxygen-air. 

2.2.  NIRS  for  measuring  AfHbOJ 

A  homodyne  frequency-domain  system  (NIM,  Philadelphia,  PA)  was  used  to  monitor  the  global 
change  of  deoxy-  and  oxy-haemoglobin  concentration  (A[Hb02j)  in  the  tumour,  as  described  previously 
(Yang  et  al  1997),  though  with  minor  modifications  to  ensure  MR  compatibility  (figure  1).  Briefly,  the 
laser  light,  which  was  emitted  from  two  laser  diodes  (at  758  nm  and  785nm),  was  amplitude-modulated  at 
140  MHz  and  time  gated  on  and  off  sequentially.  The  two  time-shared  laser  beams  illuminated  the  tumour 
surface  alternately  through  a  light-delivery  fibre  bundle  (source  fibre)  with  a  7-meter  length  and  a  3-mm 
bundle  diameter.  The  long  fibre  bundle  ensured  the  separation  of  the  NIRS  hardware  from  the  magnet. 
After  being  absorbed  and  scattered  in  the  tumour  tissue,  the  transmitted  light  was  collected  on  the 
opposite  side  of  the  tumour  by  another  fibre  bundle  of  same  length  and  diameter  (detector  fibre)  as  those 
of  source  fibre.  The  collected  optical  signal  was  then  detected  and  amplified  by  a  photomultiplier  tube 
(PMT).  The  source  and  detector  fibre  bundles  are  shown  in  figure  1 ,  and  both  of  the  probe  sheaths  at  the 
fibre  tips  were  made  of  nylon  for  MR  compatibility.  An  In-phase  and  Quadrature  (I/Q)  demodulator  chip 
was  used  to  demodulate  the  amplitude-modulated  signal  from  the  PMT. 
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In  principle,  since  the  I/Q  system  could  give  both  phase  and  amplitude  values,  we  should  be  able  to 
obtain  absolute  calculations  of  Hb02,  Hb,  and  S02  (Yang  et  al,  1997).  However,  given  the  tumor’s  small 
size  and  large  spatial  heterogeneity,  it  is  very  difficult  to  obtain  such  absolute  quantification  accurately 
using  the  conventional  algorithm  (Fishkin  and  Gratton  1993)  based  on  the  diffusion  approximation. 
Instead,  based  on  modified  Beer-Lambert’s  law,  we  can  use  the  amplitude  of  the  light  transmitted  through 
the  tumor  to  calculate  concentration  changes  in  Hb02  Hb,  and  Hb,  of  the  tumor  caused  by  respiratory 
intervention.  Namely,  changes  of  oxy-  and  deoxy-haemoglobin  concentration,  A[Hb02]  and  A[Hb], 
respectively,  can  be  derived  from  the  measured  amplitudes  at  the  two  wavelengths  (758  nm  and  785  nm) 
(Kim  et  al  2003b) 


A  [Hb02]  = 


-10.63  log(  Ah  )758  +  14.97  log(  Ah  )785 

_ A _ A 

DPF  ■  d 


A  [Hb]  = 


8.95  log(  Ah  )758  -  6.73  log(  Ab  )785 

_ A _ A 

DPF  d 


(1) 

(2) 


where  Ah  is  the  baseline  amplitude,  A,  is  the  transient  amplitude  during  measurement,  and  d  is  the  direct 
source-detector  separation.  DPF  (differential  path-length  factor)  is  a  tissue-dependent  parameter  and  is 
defined  as  the  ratio  between  the  optical  path  length  and  the  physical  separation  between  the  source  and 
detector.  Notice  that  the  units  for  A[Hb02],  A[Hb ],  and  A[Hb],otai  are  mM.  Since  DPF  is  a  variable, 
depending  on  tissue  types  and  wavelengths,  it  is  currently  difficult  to  quantify  DPF  for  tumours.  Because 
our  study  focuses  on  dynamic  changes  of  [Hb02],  A[Hb02]  values  may  be  scaled  by  a  factor  of  DPF,  in 
the  unit  of  mM/DPF,  to  obtain  characteristic  features  of  tumour  oxygen  dynamics  (Liu  et  al  2000). 


2.3.  Mathematical  model  for  blood  oxygenation  dynamics  of  tumours 

Based  on  our  previous  experimental  study,  we  derived  a  simple  mathematical  model  (Liu  et  al  2000)  to 
examine  tumour  vascular  dynamics  during  oxygen  intervention,  by  analogy  to  the  method  used  to 
quantify  regional  cerebral  blood  flow  (rCBF)  with  diffusible  radiotracers,  originally  developed  by  Kefy 
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(1951).  Accordingly,  changes  of  oxy-haemoglobin  concentration  in  tumour  vasculature  induced  by 
hyperoxic  gas  intervention,  A[Hb02],  could  be  expressed  as 

A[Hb02(t)]va™h't're  =^o['-exp(-^)]  =  41-exp(--)]  (3) 

Y  T 

where  A[Hb02]  corresponds  to  the  changes  in  oxy-haemoglobin  concentration  from  tumour  vasculature 
measured  by  the  NIRS,  Hn  is  the  arterial  oxygenation  input  for  AHb02artery  after  time  0, /is  the  blood 
perfusion  rate,  and  y  is  the  ratio  of  [Hb02]  changes  in  the  vascular  bed  to  that  in  the  veins  and  defined  as 
vascular  coefficient  of  the  tumour.  Specifically,  y  =  A[Hb02]vasculalure/A[Hb02]'e"1,  r  is  the  time  constant 
(=  y/f),  and  A=  y  H0. 

Given  the  fact  that  solid  tumours  often  develop  hypoxic  regions  which  are  poorly  perfused  in  the 
centre  as  they  grow  (Mazurchuk  et  al  1999,  Song  et  al  2002),  we  hypothesized  that  the  observed,  bi- 
phasic  feature  of  A[Hb02]  came  from  two  different  perfusion  regions  in  tumours  (well-perfused  and 
poorly  perfused  region).  Therefore,  it  is  reasonable  to  define  two  different  blood  perfusion  rates  (/i,/) 
with  two  different  vascular  coefficients  (yi,  y2)  in  the  mathematical  model.  Consequently,  eq.3  has  been 
modified  to  a  summation  of  two  exponential  expressions,  representing  two  perfusion  regions,  as 

A [Hb02 (0rc",o"'re  =  T,//0t  1  -  exp(-— )]  +  y2H0[\  -  exp(-^)] 

Y\  Y  2 

=  zf,[l-exp(-—  )]+/f2[l-exp(-— )]  (4) 

r,  r2 

where/  and  yi  are  the  blood  perfusion  rate  and  the  vasculature  coefficient,  respectively,  in  the  well- 
perfused  region,  and/  and  y2  represent  the  same  respective  meanings  in  the  poorly  perfused  region.  Also, 
it  follows  that  zf|=yi//o,  A2=  y2//0,  fi  =  yi//i  and  r2=  y2//.  Since  A],  A2,  r,,  and  r2can  be  determined  by 
curve-fitting  equation  (4)  to  the  dynamic  NIRS  measurements,  we  can  obtain  the  ratios  for  two  vascular 
coefficients  and  two  blood  perfusion  rates,  as  follows  (Liu  et  al  2000): 
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ZL=Aj_  fj  A|/r|  ^ 

Y2  A2  ’  f2  A  2  / tj_ 

These  two  ratios  enable  us  to  understand  more  about  tumour  vascular  structures  and  blood  perfusion  rates. 
Namely,  yxly2  may  be  associated  with  the  vascular  volume  fraction  between  the  two  regions,  and  flf2 
reflects  the  perfusion  ratio  between  the  two  regions.  For  example,  if  flf  is  substantially  greater  than  1, 
the  two  compartments  have  significantly  different  perfusion  rates. 

This  mathematical  model  has  been  utilized  well  to  explain  our  experimental  animal  studies  (Liu  et  at 
2000,  Gu  et  a!  2003,  Kim  et  al  2003b).  Moreover,  it  has  been  supported  by  computer  simulations  (Kim 
and  Liu  2005)  and  phantom  study  (Kim  and  Liu  2004)  that  the  bi-phasic  or  bi-exponential  feature  can  be 
well  present  if  both  a  slow  and  a  fast  perfusion  of  flow  co-exist  within  the  interrogated  area  or  volume  by 
the  N1R  source  and  detector.  In  this  paper,  our  theoretical  analysis  is  also  based  on  the  same  mathematical 
model,  as  given  in  equations  (4)  and  (5). 

2. 4.  FREDOM  for  measuring  pOi 

MRI  experiments  were  performed  using  a  Varian  Inova  4.7  T  horizontal  bore  system  equipped  with 
actively  shielded  gradients.  Shimming  was  performed  on  the  tumour  tissue  water  signal  to  reduce  the 
line-width  to  less  than  100  Hz.  ‘H  MRI  (200.1  MHz)  Tl-weighted  reference  images  were  acquired  with 
TR/TE  of  150/10  ms  and  40x40  mm2  field  of  view.  Following  'H  MRI,  corresponding  19F  MR  images 
(188.3  MHz)  with  a  matrix  size  of  32x32  (1.25  mm  per  pixel)  were  obtained  to  show  the  distribution  of 
HFB  in  the  tumour.  A  single  4-mm  thick  slice  was  chosen  to  include  all  of  the  injected  HFB  as  the  plane 
of  injection  may  not  be  perfectly  aligned  with  the  imaging  plane.  The  FREDOM  (Fluorocarbon 
Felaxometry  by  £cho-planar  imaging  for  Dynamic  Oxygen  Mapping)  approach  was  used  to  measure  p02, 
as  described  in  detail  previously  (Hunjan  et  al  2001).  The  spin  lattice  relaxation  rate  (Rl=  1/T1)  of  HFB 
is  highly  sensitive  to  changes  in  p02,  but  not  sensitive  to  temperature  variations:  a  deviation  of  1  °C  in 
temperature  will  introduce  a  deviation  of  only  0.13  torr  in  p02  estimate,  when  p02  is  about  5  torr  (Zhao  et 
al  2004).  T1  maps  were  computed  on  a  voxel-by-voxel  basis  using  non-linear  least-squares  data  fitting 
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by  the  Gauss-Newton  method.  We  applied  a  threshold  to  the  raw  T1  data  in  order  to  remove  random 
noise,  i.e.,  voxels  with  T1  error  >3.6  sec  or  T1  error/Tl  >50%  were  disregarded.  Maps  of  p02  values  with 
1.25  mm  pixel  size  were  obtained  from  the  T1  maps  using  the  equation,  p02  =  (l/Tl-0.0835)/0. 001876 
(Hunjan  et  al  2001). 

Each  p02  map  took  6  and  1/2  minutes  to  acquire.  Three  baseline  p02  data  sets  were  acquired  over  24 
mins  for  all  tumours,  while  the  rats  breathed  air,  after  which  the  rats  were  repeatedly  exposed  to  oxygen 
or  carbogen  (95%  02  and  5%  C02)  interventions.  Five  p02  maps  were  obtained  during  each  subsequent 
gas  switch  period.  Typically,  for  a  five  gas-intervention  sequence  {eg.,  air-oxygen-air-carbogen-air),  a 
total  of  23  p02  maps  were  obtained  over  a  period  of  3  h.  Due  to  thresholding,  some  voxels  did  not  appear 
in  all  p02  maps  and  these  were  discarded.  For  temporal  analysis,  voxels  were  selected  as  only  those 
which  provided  consistently  reliable  data  {i.e.,  satisfied  all  the  thresholding  criteria  specified  above)  for 
all  23  measurements  over  the  time  course  with  a  range  of  5  to  44  acceptable  voxels  for  the  ten  tumours. 
The  slope  of  dynamic  p02  changes  (rate)  was  defined  as  d(p02)/dt  and  d(p02’)/dt  in  response  to 
increasing  or  decreasing  inhaled  F02  (Fraction  of  02),  respectively. 

2.5.  Statistical  Analysis 

Linear  regression  analysis  was  used  to  calculate  the  correlation  between  the  NIRS-derived  tumour 
hemodynamic  parameters  {i.e.,  A[Hb02],  A\!t i,  A2It2i  f/fi)  and  the  F7?£D0M-determined  tumour 
parameters  {i.e.,  p02,  d(p02)/dt,  d(p02’)/dt).  Data  are  presented  as  mean  ±  standard  deviation  (SD)  and 
paired  Student-t  tests  compared  the  effects  of  oxygen  and  carbogen  on  A[Hb02]  and  p02. 

3.  Results 

3. 1.  Dynamic  response  of  A[Hb02]  measured  by  NIRS 

Figure  2  shows  a  group-averaged,  temporal  profile  of  A[Hb02]  from  seven  13762NF  rat  breast  tumours 
with  air-oxygen-air-carbogen-air  intervention,  displaying  apparent  biphasic  responses  to  both 
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interventions.  Both  single  (eq.  3)  and  double-exponential  (eq.  4)  curve  fitting  were  tested  for  the  carbogen 
intervention  in  one  representative  tumour  (1.6  cm3),  as  in  figure  3 A.  The  maximal  A[Hb02]  achieved 
with  oxygen  challenge  was  compared  with  that  of  carbogen,  and  revealed  no  significant  difference 
between  oxygen  and  carbogen  intervention  (p>0.3);  indeed,  there  was  a  strong  correlation  between  the 
maximal  A[Hb02]  values  with  these  two  interventions  (R2>0.75,  figure  3B),  consistent  with  our  previous 
observations  (Gu  et  al  2003).  No  correlation  was  found  between  the  perfusion  rate  ratio  (/!//>)  and  tumour 
size  (R2=0.16).  Vascular  oxygen  dynamics  in  response  to  interventions  are  provided  for  individual 
tumours  as  supplementary  material. 

3.2.  p02  measurements  by  F RE  DOM 

Overlay  of  l9F  and  'H  MR  images  demonstrated  that  HFB  was  distributed  in  both  central  and  peripheral 
regions  of  tumour  (not  shown).  Individual  p02  values  taken  from  each  voxel  ranged  from  hypoxia  (<  1 
torr)  to  35  torr  under  baseline  conditions  (figure  4).  Mean  baseline  p02,  which  are  averaged  over  all 
voxels  in  the  first  three  p02  maps,  ranged  from  hypoxia  (<  5  torr)  to  27  torr  with  a  hypoxic  fraction  (HF5; 
fractional  voxels  that  are  less  than  5  torr)  ranging  from  0  to  100%  (mean  36%)  and  summarised  in  table  1. 
A  strong  correlation  was  found  between  mean  baseline  p02  and  HF5  (R2>0.85,  figure  5).  Administration 
of  oxygen  or  carbogen  produced  significant  increases  in  tumour  p02,  as  shown  in  the  p02  maps,  graphs 
and  table.  The  tumour  p02  values  shown  here  are  volume-averaged  over  the  entire  slice.  The  p02 
responses  to  respiratory  challenge  from  the  group  of  seven  tumors,  measured  simultaneously  by 
FREDOM  and  NIRS,  are  also  shown  in  figure  2.  Baseline  measurements  (breathing  air)  were  generally 
stable,  and  altering  the  inhaled  gas  to  oxygen  or  carbogen  induced  rapid  and  significant  changes  in  both 
p02  and  A[Flb02]  (p<0.001)  Upon  return  to  air  (baseline),  A[Hb02]  dropped  quickly  and  significantly 
within  16  mins,  and  then  more  slowly,  for  the  next  24  mins,  whereas  the  p02  decrease  was  more  gradual. 
Altering  the  inhaled  gas  to  carbogen  also  produced  a  rapid  increase  in  both  p02  and  A[Hb02].  Upon 
return  to  air  breathing  from  carbogen,  both  A[Hb02]  and  p02  showed  a  similar  trend  to  that  following 
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oxygen.  As  expected,  all  ten  tumours  showed  a  significant  increase  in  p02  (volume  average  over  entire 
slice),  and  decrease  in  hypoxic  fraction  (HF)  in  response  to  oxygen  or  carbogen  inhalation.  The 
magnitude  of  response  to  either  hyperoxic  gas  was  correlated  (R2>0.79),  as  was  the  maximum  volume- 
averaged  p02  achieved  with  either  gas  (R2>0.83).  The  rate  of  increase  with  oxygen  challenge,  d(p02)/dt, 
was  significantly  faster  than  the  return  to  baseline,  d(p02’)/dt,  for  oxygen  intervention  (p<0.02),  but  no 
difference  was  observed  for  carbogen  (p>0.1,  figure  6).  The  mean  p02  values  averaged  over  the  final 
three  p02  maps  during  exposure  to  carbogen  breathing  were  significantly  higher  than  oxygen  (p<0.01), 
and  the  tumour  hypoxic  fraction  was  generally  eliminated  during  carbogen  breathing  (n  =  7  of  9  tumours, 
table  1). 

3.3.  The  relationship  between  p02  and  A[Hb02]  with  respect  to  hyperoxic  gas 

Taken  as  a  group  of  10  tumours,  there  was  no  apparent  relationship  between  the  magnitude  of  the  change 
in  tumour  vascular  oxygenation  (A[Hb02])  and  change  in  p02  (R2<0.1).  However,  if  tumours  were 
divided  into  two  sub-populations,  then  two  separate  correlations  were  found  each  with  similar  slope 
(figure  7A).  There  was  also  a  correlation  (R2>0.7)  between  the  perfusion  rate  ratio,  ft  fa,  derived  from 
fitting  the  A[Hb02]  curve  and  the  mean  p02  values  achieved  with  hyperoxic  gas  intervention  (figure  7B). 
Assessment  of f/fi  is  predicated  on  biphasic  behaviour  with  respect  to  interventions,  which  was  observed 
in  most  cases  (13  of  16  measurements).  There  was  also  a  positive  correlation  between  A\/t\  (the  fast 
component  of  biphasic  A[Hb02])  and  the  d(p02)/dt  of  well  oxygenated  voxels  ( i.e .,  those  with  p02  values 
>  10  torr  under  oxygen  or  carbogen  intervention)  (R2>0.5,  figure  8).  However,  no  correlation  was  found 
between  d(p02)/dt  and  A2/t2  (the  slow  component). 

4.  Discussion 

Integration  of  diverse  imaging  techniques  can  be  technically  challenging,  since  each  modality  has  specific 
technical  constraints  and  requirements.  Here,  the  NIR  system  was  modified  to  use  longer  optical  fibres 
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and  any  metal  components  were  eliminated  from  the  fibre  tips.  Due  to  the  spatial  restrictions  within  the 
bore  of  the  magnet  a  Helmholtz  coil  was  built  specifically  providing  access  to  both  the  tumour  and  fibres. 
The  fibres  required  sufficient  flexibility  to  allow  them  to  be  bent  through  requisite  angles  within  the 
confines  of  the  magnet  bore.  To  date  there  have  been  few  reports  of  simultaneous  data  acquisition  by  NIR 
and  MRI  on  phantom  (Pogue  et  al  2003),  human  brain  (Toronov  et  al  2001 ,  Chen  et  a!  2003)  and  human 
breast  (Ntziachristos  et  al  2000,  Gulsen  et  al  2002,  Ntziachristos  et  al  2002,  Gu  et  al  2004a). 

In  the  present  study,  global  average  A[Hb02]  was  measured  by  NIRS,  and  p02  maps  were  obtained 
simultaneously  by  l9F  MRI.  We  used  transmission  mode  NIRS  in  order  to  interrogate  deep  tumour  tissue. 
While  NIRS  is  a  global  measurement,  the  region  sampled  by  it  is  predominantly  a  banana-shape  region 
joining  the  locations  of  the  source  and  detector  (Arridge  1995,  Arridge  and  Schweiger  1995),  and  we 
recently  showed  that  typically  15  to  30%  of  the  vascular  volumes  of  rat  tumours  are  interrogated  by  NIR 
using  this  configuration  (Gu  et  al  2004b).  We  positioned  the  NIRS  probes  (source  and  detector)  to  be  in 
the  same  plane  as  the  HFB  injection  to  ensure  maximum  overlap  between  the  regions  sampled  by  the  two 
techniques.  Utilising  our  previously  developed  mathematical  model  (Liu  et  al  2000),  multiple  hemo¬ 
dynamic  parameters  were  derived  for  A[Hb02]  (A  d T\A2 / r2  and/]//))  to  be  compared  with  p02.  Our  results 
demonstrate  that  oxygenation  parameters  measured  from  both  techniques  show  significant  and  consistent 
elevation  in  tumour  oxygenation  during  the  hyperoxic  gas  interventions.  As  reported  previously,  the 
magnitude  of  the  vascular  response  was  similar  with  both  hyperoxic  gases  (Gu  et  al  2003).  As  expected, 
A[Hb02]  increased  much  faster  than  p02  in  all  ten  tumours,  indicating  that  change  in  tumour  vascular 
oxygenation  precedes  tumour  tissue  oxygenation.  This  observation  is  consistent  with  our  previous  studies 
in  this  tumour  type  measured  simultaneously  by  NIRS  and  fibre-optic  probes  (Gu  et  al  2003),  as  well  in 
the  Dunning  prostate  R3327-AT1  tumours  measured  sequentially  by  NIRS  and  l9F  MRI  (Kim  et  al 
2003b). 

We  have  previously  demonstrated  the  application  of  FREDOM  to  monitor  tumour  oxygen  dynamics 
in  diverse  rat  prostate  tumours  (Zhao  et  al  200 1 ,  Zhao  et  al  2002),  human  tumour  xenografts  (Mason  et 
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al  2002)  and  a  few  breast  tumours  (Song  et  al  2002).  Here,  mean  baseline  was  p02=l  2  ±  1 0  torr  for  the  1 0 
tumours,  which  is  lower  than  reported  previously  (Song  et  al  2002),  but  entirely  consistent  with  the  newer 
anaesthetic  protocol  (air  or  21%  oxygen,  as  opposed  to  33%  02  previously).  There  was  a  strong 
correlation  between  baseline  p02  and  hypoxic  fraction  (figure  5),  as  we  have  previously  found  using 
Dunning  prostate  R3327-HI  tumours  (Zhao  et  al  2001).  The  p02  achieved  with  carbogen  in  this  study  was 
significantly  higher  than  with  oxygen  and  carbogen  appeared  to  be  more  effective  at  eliminating  the 
hypoxic  fraction.  However,  carbogen  was  generally  applied  second  in  our  experimental  protocols,  and  it 
is  highly  likely  that  the  initial  oxygen  primed  the  tumour.  Indeed,  while  oxyhaemoglobin  generally 
returned  to  baseline  during  the  air  breathing  episode  between  hyperoxic  gases,  it  is  clear  that  p02 
remained  elevated  (figure  2).  Both  Ap02  and  the  maximum  p02  achieved  with  either  gas  were  closely 
correlated.  Similar  behaviour  was  reported  previously  based  on  measurements  using  fibre  optic  probes 
(Gu  et  a I  2003).  For  the  group  of  tumours  in  this  study,  mean  baseline  p02  did  not  provide  a  good 
indication  of  response  to  hyperoxic  gas  (Ap02  or  p02max).  However,  considering  the  fate  of  individual 
voxels  in  separate  tumours  with  respect  to  intervention,  strong  linear  correlation  was  observed  in  some 
tumours  (tumour  #3,  4,  6  and  9)  between  initial  mean  baseline  p02  and  the  maximum  p02  at  the  same 
location  (voxel)  in  response  to  carbogen  or  oxygen  breathing  consistent  with  the  results  of  Song  et  al 
(2002).  Other  tumours  showed  no  correlation. 

The  rate  of  p02  response  to  either  gas  was  similar  (figure  6),  but  decrease  upon  return  to  air  was 
significantly  faster  in  the  case  of  oxygen,  but  not  carbogen.  While  the  relationships  between  Ap02  and 
A[Hb02]  are  not  obvious,  there  exists  a  significant  correlation  between  flf2  and  mean  p02  values  achieved 
with  hyperoxic  gas  intervention.  Most  flf2  values  are  between  5  and  60  (figure  7B),  implying  that  the 
blood  perfusion  rate  in  the  well  perfused  region  is  much  higher  than  that  from  the  poorly  perfused  region. 
Not  surprisingly  the  higher  the  ratio  of  well-perfused  to  poorly-perfused  regions,  the  higher  the  mean  p02 
values  achieved  by  the  hyperoxic  gas  interventions. 
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Previous  studies  have  demonstrated  that  tumour  tissue  oxygenation  could  be  strongly  affected  by 
changes  in  tumour  blood  flow  measured  locally  by  laser  Doppler  flowmetry  (Vaupel  et  al  1994).  It  was 
also  reported  that  spontaneous  fluctuations  in  flow  and  peri-vascular  p02  are  correlated  at  the  on  micro- 
regional  level  (Kimura  et  al  1996,  Braun  et  at  1999).  In  our  study,  we  are  able  to  derive  the  perfusion  rate 
ratios  between  the  well  perfused  and  poorly  perfused  regions  from  the  haemoglobin  concentration 
responses  to  the  intervention  measured  by  NIRS.  We  found  that  there  was  a  significant  correlation 
between  d(p02)/dt  and  A \/t\  (f\  proportional  to  A but  not  d(p02)/dt  and  Ad r2  (/)  proportional  to  Ad r2), 
provided  that  the  p02  readings  were  selected  from  well  oxygenated  or  responsive  voxels  (figure  8).  This 
linear  correlation  suggests  that  the  rate  of  change  in  p02  is  closely  related  to  the  perfusion  rate  in  the  well 
perfused  region,/].  In  other  words,  the  dynamic  changes  in  p02  of  those  regions  responsive  to  hyperoxic 
gas  intervention  may  be  attributed  to  fast  tumour  vascular  perfusion,  rather  than  to  the  slow  perfusion  in 
tumour  vasculature.  We  believe  these  results  provide  a  valuable  association  between  tumour  vascular 
oxygenation  and  tumour  p02  determined  simultaneously  by  the  optical  and  NMR  measurements. 

A  goal  had  been  to  develop  a  low-cost,  simple,  fast  surrogate  measurement  of  p02  based  on  NIRS  of 
the  oxygenation  status  of  endogenous  haemoglobin.  Both  figures  7B  and  8A  suggest  that  there  exist  linear 
relationships  of  p02  with  the  NIRS  measurable  parameters  fi/f2  and  with  A/ti  that  are  measurable  for 
NIRS.  However,  the  correlation  between  Ap02  and  A[Hb02]  (figure  7A)  seems  to  be  more  complex  with 
a  separation  of  the  tumours  into  two  groups.  It  may  be  noteworthy  that  the  majority  (6  of  7)  of  the 
tumours  associated  with  the  correlation  on  the  left  hand  side  of  the  graph  had  a  high  initial  global  p02 
(i.e.,  mean  baseline  p02  >5  torr),  whereas  2  out  of  3  of  the  tumours  on  the  right  had  low  initial  global  p02. 
A  possible  interpretation  relates  to  change  in  the  shape  of  the  haemoglobin-oxygen  dissociation  curve  (the 
Hill  curve).  A  given  p02  response  can  correspond  to  different  changes  in  haemoglobin  saturation 
depending  on  where  the  change  occurs  on  the  Hill  curve.  At  lower  initial  p02  there  may  be  a  substantial 
Ap02  with  little  A[Hb02].  By  contrast,  at  higher  initial  p02  the  same  increase  in  p02  will  produce  a 
greater  effect  on  the  saturation.  While  A[Hb02]  is  proportional  to  the  change  of  S02  (with  the  assumption 
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of  little  change  in  total  haemoglobin  concentration  during  intervention),  two  subpopulations  of  tumours  in 
the  relationship  between  Ap02  and  A[Hb02]  could  be  the  result  of  different  initial  p02.  Similarly,  two 
subpopulations  of  tumours  were  observed  in  the  relationship  between  initial  S02  and  the  carbogen- 
induced  change  in  saturation  in  the  study  of  Hull  et  al  (1999).  Such  an  effect  also  confounds  the  direct 
correlation  of  BOLD  MRI  response  to  changes  in  p02  (Baudelet  and  Gallez  2002).  Significantly, 
preliminary  data  reported  by  Gu  et  al  (2004a)  for  simultaneous  BOLD  MRI  and  NIR  in  tumors  showed  a 
strong  linear  correlation  in  response  to  hyperoxic  gas  challenge,  as  also  reported  by  chen  et  al  (2003)  in 
the  rat  brain. 

A  major  concern  is  tumour  heterogeneity,  as  recognised  throughout  the  literature  and  shown  here  by 
l9F  MRI  (figure  4).  Indeed,  we  have  obtained  some  preliminary  data  using  a  single  NIR  source  and  three 
detectors  placed  on  various  regions  across  a  tumour  (Kim  et  al  2003a),  showing  that  each  region  of  the 
tumour  responded  differently  to  hyperoxic  gas,  in  terms  of  the  extent  and  rate,  indicating  the 
heterogeneity  of  tumour  vasculature.  Spatial  discrimination  will  be  even  more  critical,  if  such  studies  are 
transferred  to  human  breast  cancer,  where  the  tumour  is  surrounded  by  normal  tissue  (Brooksby  et  al 
2003).  Nevertheless,  we  believe  this  haemodynamic  model  and  correlation  between  tumour  vascular 
oxygenation  and  p02  provides  valuable  insight  into  the  tumour  compartment  of  such  a  mixed  system  and 
explores  dynamic  signatures  of  breast  tumours,  which  could,  in  turn,  enhance/assist  human  breast  cancer 
diagnosis  and  prognosis. 

In  summary,  by  studying  tumour  vascular  oxygenation  concomitantly  with  changes  in  tumour  oxygen 
tension,  we  found  several  significant  correlations  between  rates  and  magnitudes  of  vascular  and  tissue 
responses.  This  study  also  demonstrates  the  feasibility  of  conducting  simultaneous  NIRS  and  MRI 
oximetry.  We  believe  the  correlation  of  tumour  vascular  oxygenation  and  tumour  tissue  p02  can  provide 
valuable  insights  into  tumour  pathophysiology  and  response  to  interventions. 
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Figure  Legends 
Figure  1 

Schematic  of  experimental  set-up.  Z  is  along  the  bore  of  the  magnet,  and  X  along  the  axis  of  the  rf  coil. 
PMT  represents  a  photomultiplier  tube,  I/Q  represents  an  In-phase  and  Quadrature  (I/Q)  demodulator  chip 
and  ADC  represents  an  analogue  to  digital  converter. 

Figure  2 

Temporal  profiles  of  A[Hb02]  (curve  only)  and  mean  p02  (curve  with  ♦)  in  response  to  respiratory 
challenge.  Both  A[Hb02]  and  p02  are  group-averaged  data  from  seven  13762NF  rat  breast  tumours  with 
an  intervention  sequence  of  air-oxygen-air-carbogen-air,  measured  simultaneously  by  NIRS  and 
FREDOM.  The  standard  deviations  of  the  volume-averaged  baseline  p02  from  7  tumours  are  indicated  in 
figure  2.  The  unit  for  A[Hb02]  is  mM/DPF. 

Figure  3 

A)  Representative  dynamic  responses  of  A[Hb02]  to  gas  intervention  (carbogen).  The  unit  for 
A[Hb02]  is  mM/DPF.  Single  exponential  curve  fitting  yielded  A[Hb02]  =  0.22{1  -  exp[-(t- 

106.8) /2.5]}  (R2=0.64),  and  double  exponential  fitting  resulted  in  A[Hb02]  =  0. 15 { 1  -  exp[-(t- 

106.8) /0.58]}+  0.13{1  -  exp[-(t-106.8)/23.2]}  (R2=0.88). 

B)  Relationship  of  maximum  A[Hb02]  (in  mM/DPF)  in  breast  tumours  in  response  to  switching 
from  air  to  oxygen  and  to  carbogen  (R2>0.75). 


Figure  4 

Representative  p02  maps  (1.25  mm  resolution)  obtained  using  FREDOM  overlaid  on  the  !H  anatomic 
image  (FOV=  4  cm  x  4  cm).  Left)  Rat  breathing  air.  Centre)  breathing  100%  02  (fifth  map  after 
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switching  from  air),  Right)  breathing  carbogen  (95%  02  +  5%  C02)  (fifth  map  after  switching  from  air). 
The  p02  maps  show  distinct  heterogeneity. 

Figure  5 

Correlation  between  baseline  p02and  hypoxic  fraction  (HFs)  measured  using  FREDOM  (R2>0.85). 

Figure  6 

Mean  d(p02)/dt  (open)  and  d(p02’)/dt  (shaded)  j^SD  for  eight  tumours  with  both  interventions,  when  gas 
was  switched  from  air  to  the  hyperoxic  gas  and  back  to  air,  respectively.  p02  is  the  mean  value  of  all 
acceptable  voxels  appearing  in  the  five  maps  during  the  oxygen  or  carbogen  intervention.  d(p02)/dt  is  the 
slope  of  regression  line  of  five  p02  readings  versus  time  during  oxygen  or  carbogen  interventions,  and 
d(p02’)/dt  is  the  slope  of  regression  line  for  five  p02  readings  versus  time  during  air  interventions  after 
oxygen  or  carbogen.  The  rates  showed  a  significant  difference  with  oxygen  (p<0.02),  but  not  with 
carbogen  (p>0.1). 

Figure  7 

A)  Correlation  between  maximum  A[Hb02]  and  change  in  p02  with  respect  to  hyperoxic  gas  intervention 
for  two  groups  of  tumours  (•)  (group  1 :  #1,  2,  6,  7,  8,  9,  10;  R2  >  0.51)  and  (A)  (group  2:  #3,  4,  5;  R2 
>0.82).  The  unit  for  A[Hb02]  is  mM/DPF. 

B)  Correlation  between  mean  p02  achieved  with  hyperoxic  gas  breathing  and  perfusion  rate  ratio  (/!//)) 
for  tumours  with  biphasic  response  to  intervention  (R2>0.7).  p02  is  the  mean  value  for  the  final  three 
p02  maps  under  hyperoxic  intervention,  selected  from  all  voxels  appearing  in  the  final  three  p02 
maps  during  oxygen  (♦)  or  carbogen  (□)  intervention. 

Figure  8 

d(p02)/dt  vs  A  |/ri  determined  from  A[Hb02]  for  tumours  with  biphasic  response  to  interventions,  showing 
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a  positive  correlation  (R2>0.5).  p02  are  obtained  from  the  mean  value  of  all  well  oxygenated  voxels  (i.e., 
the  maximum  p02  >  10  torr  under  oxygen  or  carbogen  intervention)  appearing  in  the  five  p02  maps 
during  oxygen  (♦)  or  carbogen  (□)  intervention 
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Figure  1 
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Figure  3A 
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Table  1.  Oxygen  Tension  (p02)  in  ten  rat  mammary  13762NF  adenocarcinomas.  p02  in  baseline  is  the  average  value  of  all  voxels 
in  the  first  three  maps,  while  p02  for  oxygen  or  carbogen  is  the  average  value  of  all  voxels  in  the  final  three  p02  maps  during 
exposure  to  the  hyperoxic  gas.  Hypoxic  fraction  (HF)  is  the  percentage  of  voxels  with  p02  values  less  than  5  torr  (HF5)  or  10  torr 
(HFio)  and  HF5  and  HF|0  are  the  average  value  of  final  three  maps  in  baseline,  oxygen  or  carbogen  respectively.  The  mean  p02 
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Table  2.  Summary  of  vascular  oxygen  dynamics  for  the  experimental  tumours.  Two  amplitudes 
(A  |,  A2 )  and  two  time  constants  (ru  r2)  were  determined  by  curve-fitting  the  dynamic  NIRS 
measurements  using  a  double-exponential  expression.  Nine  out  of  10  tumours  were  observed  to 
have  double-exponential  features  with  either  oxygen  or  carbogen  intervention.  A\  is  significantly 
smaller  than  A2  with  oxygen  (p<0.01),  but  no  significant  differences  in  carbogen  (p>0.19). 


Tumour _ Double  exponential  fitting 


No. 

intervention 

A\  (AIJ) 

^2(AU) 

T\  (min) 

r2  (min) 

R 

a1 

oxygen 

0.11  ±0.002 

0.12  ±0.002 

0.49  ±0.02 

6.6  ±0.2 

carbogen 

0.15  ±0.001 

0.13  ±0.002 

0.58  ±  0.01 

23.2 ±  1.1 

0.94 

a2 

oxygen 

0.03  ±0.002 

0.1 1  ±0.005 

0.46  ±0.18 

25.6±  2.9 

0.91 

carbogen 

— 

— 

— 

— 

— 

a3 

oxygen 

0.15  ±0.001 

0.38  ±0.003 

0.31  ±0.02 

20.1  ±0.39 

carbogen 

0.23  ±0.005 

0.2  ±0.004 

1.2  ±0.04 

10.8  ±0.45 

0.95 

“4 

oxygen 

— — 

— 

— 

— 

— 

carbogen 

0.16  ±0.002 

2.3  ±0.27 

0.29  ±0.02 

155.8  ±  21 

a5 

oxygen 

0.12  ±0.002 

0.23  ±  0.002 

0.55  ±0.03 

11. 6  ±0.27 

0.94 

carbogen 

— 

— 

— - 

— 

— — 

% 

oxygen 

0.06  ±0.001 

0.26  ±0.005 

0.13  ±0.01 

31.2  ±  1.1 

0.96 

carbogen 

0.03  ±0.001 

0.08  ±0.001 

0.04  ±0.02 

13.9  ±  0.6 

0.78 

a? 

oxygen 

— 

— — 

— 

~~ 

— — 

carbogen 

b8 

oxygen 

0.11  ±0.002 

0.2  ±0.001 

0.24  ±0.02 

7.1  ±0.1 

carbogen 

0.2  ±0.002 

0.12  ±0.001 

0.49  ±0.01 

9.1  ±0.2 

c9 

carbogen 

0.01  ±0.0004 

0.02  ±0.0004 

0.14  ±0.02 

3.9±  0.12 

0.64 

d10 

oxygen 

0.05  ±0.002 

0.03  ±  0.002 

0.46  ±0.05 

5.6  ±0.56 

0.68 

Mean 

±SD 

.... 

0.11  ±0.068 

0.32  ±0.59 

0.42  ±0.3 

24.9  ±40.2 

.... 

a:  air->02  ~>air  ->carbogen~>air. 
b:  air  ->  carbogen  ->  air  02  ->  air. 
c:  air  ->  carbogen  air. 
d:  air  oxygen  ->  air. 
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